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In this paper, we investigate the Ricci dark energy model with perturbations through the joint 
constraints of current cosmological data sets from dynamical and geometrical perspectives. We use 
y—^ • the full cosmic microwave background information from WMAP seven-year data, the baryon acoustic 

^— ( ' oscillations from the Sloan Digital Sky Survey and the Two Degree Galaxy Redshift Survey, and 

, type la supernovae from the Union2 compilation of the Supernova Cosmology Project Collaboration. 

^Nj . A global constraint is performed by employing the Markov chain Monte Carlo method. With 

' the best-fitting results, we show the differences of cosmic microwave background power spectra 

O , and background evolutions for the cosmological constant model and Ricci dark energy model with 

• perturbations. 
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I. INTRODUCTION 



u: 

The present observations from WMAP seven-year cosmic microwave background (CMB) radiation measurements 
^p-i large-scale structure surveys 0, 111, and updated type la supernovae (SNIa) have confirmed at a high confidence 

Q , level the first implication from SNIa in 1998 01 that our Universe is undergoing an accelerated expansion. It has 
i—{ ■ been an important issue of modern cosmology to unveil the mysterious face of the driving force of current accelerated 
c/3 , expansion. One of the theoretical viewpoints is that dark energy with negative pressure drives current expansion. A 
flood of dark energy models have been proposed, such as a natural cosmological constant ( ACDM model) , dynamical 
dark energy models with scalar fields ^ or with exotic equation of state [1, [l3l , and so on. Please see hj for an 
up-to-date review on dark energy. 

Although the cosmological constant model keeps a good fit to the current observations, it suffers from the fine- 
tuning problem and the cosmic coincidence problem. In the final analysis, both are related to the energy density of 
vacuum energy. Therefore, the cosmological constant problem is in essence an issue of quantum gravity. Before the 
complete theory of quantum gravity is established, the dark energy model on the basis of the holographic principle 
of quantum gravity theory was built in Ref. [loj by applying the energy bound proposed by Cohen, Kaplan, and 
Nelson [l^j to cosmology. Namely, for a system with size L and UV cutoff A without decaying into a black hole, it is 
required that the total energy in a region of size L should not exceed the mass of a black hole of the same size; thus 
L^PA < LMpi. It reveals a duality between the UV cutoff and IR cutoff. In cosmology, the UV cutoff is related to 
the dark energy density, and the IR cutoff is related to the large scale of the Universe. The holographic dark energy 
models with different IR cutoffs have been studied and the detailed discussions can be seen in Refs. fl3,[lMl. 
In Ref. ^13], Gao, Chen, and Shen first proposed that dark energy density is in direct proportion to the Ricci 
■ scalar curvature, R, and referred to this as the Ricci dark energy (RDE) model. The distinct characteristic of the 
RDE model is that the cosmic coincidence problem and causality problem can be naturally solved. However, the 
physical motivation for such a model was vague in Ref. [l^. SulDScquently, Cai, Hu, and Zhang J^T)] were inspired 
by the motivation of the causal entropy bound in the framework of holography, found the causal connection scale, 
RcCi and put forward an appropriate physical mechanism for this model. In [l7|, it was pointed out that the scale 
on which the black hole can be formed must be within Rcc, which is set by the "Jeans" length of the perturbations 
to i?^^ = Max{H + 2H'^, -H), while one assumes a black hole in the Universe is formed by gravitational collapse 
of perturbations of cosmological spacetime. Only if one takes the causal connection scale with R^^ = H + 2H^ as 
the IR cutoff, can the dark energy model be consistent with current observations [l3l- It is found that Rqq oc R. 
Therefore, the RDE is stressed to be the holographic dark energy with the IR cutoff of the causal connection scale. 
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The RDE density reads 

Pde = SM^ia{H + 2H'] , (1) 



where Mpi is the reduced Planck mass Mpi = 1/ VSvrG, and a is the dimensionless parameter and has been constrained 
by using the distance measurements from different joint observations 18, 19]. In particular, in these studies, the 
compressed CMB information is used, including the "shift parameter," the "acoustic scale," and the photon decoupling 
epoch, the values of which are obtained on the basis of a given model in advance [1^. That is to say, before these 
distance parameters are used to constrain a model, one should renew to get their values basing on the constrained 
model. In this paper, we revisit the cosmological parameters in the RDE model by using the full CMB power spectrum 
data from WMAP7 and additional distance measurements from baryon acoustic oscillations (BAOs) and Union2 SNIa 
including 557 data points. Compared with the previous works, we not only focus on the constraint of background 
evolution, but also take dark energy perturbations into account. For there will be dark energy perturbations according 
to the perturbed conservation equations when the dark energy is not a cosmological constant. What is more, the 
perturbations of dark energy play an important role when the dark energy model is confronted with the current 
observations jllj . The inclusion of dark energy perturbations leaves an imprint on the CMB power spectrum ^23, [23j 
and leads to the changes of parameters' spaces in the complete fitting to observational data [23-28]. In particular, 
there are the significant changes on the large-scale CMB power spectra by the late integrated Sachs- Wolfe (ISW) 
effect [29| . which arises from a time-dependent gravitational potential when dark energy dominates in the Universe 
at late time for a flat dark energy model. The magnitude of the ISW-power spectrum is dependent on dark energy 
models and dark energy perturbations [13, H^, l3Cll - l35j ]. Therefore, a complete constraint on the RDE model with 
perturbations is performed by employing the Markov chain Monte Carlo (MCMC) method. 

The paper is organized as follows. In the next section, we review the background equations for the RDE model and 
the perturbation equations in the conformal Newtonian gauge. In Sec. Ill, we perform a global fitting to observational 
data by using the MCMC method and discuss the constraint results. The last section is the conclusion. 

II. REVIEW OF THE DYNAMICS OF RICCI DARK ENERGY: BACKGROUND AND 

PERTURBATION EVOLUTION 

First of all, we give a brief review on the background evolution of the RDE model in the framework of the unper- 
turbed and spatially flat Friedmann-Robertson- Walker metric. The Einstein field equation is written as 



H^^^{Pm+Pde): (2) 



where the energy density of RDE, pde, is rewritten with the variable x = Ina as 



' 1 dH"^ 

Pde = 3aAf^, ( -— + 2i/2) . (3) 



After using the dimensionless definitions E = jf- and i},no = ta^tttt?, we can rewrite Eq. (151) as 



E'^il^oe-'^ + a(^^^+2E^y (4) 
The solution of this first order differential equation with the initial condition E^{x — 0) = 1 is obtained: 

Z — a \ 1 — a J 

From Eqs. pisp . we can get the dimensionless energy density of RDE 



2- a \ 2- a 
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where fldeo is the present value of the dimensionless energy density of RDE. Then the corresponding equation of state 
(EOS) of RDE is obtained by the conservation equation 



Idlnfldeix) 
w[x) = — 1 — - 



3 dx 



(7) 



When the EOS of dark energy is not constantly equal to -1, it is necessary to consider dark energy perturbations. 
We work in the conformal Newtonian gauge ^, the perturbed metric being 

ds^ = -a(?7)2{(l + 2i:)dr]^ - (1 - 2(?!))%dxMx^}. (8) 

The perturbed conservation equations of dark energy perturbations [S^ read 

6' = -{l + w){e-34>')-3n{cl~w)5, (9) 

/ 9 
7/7 C 

e' = -n{i-3w)9 e + —^k^6 + k^^, (lo) 

1 + w 1 + w 

where 6 is the energy density perturbation, and 9 is the velocity divergence perturbation and is related to the velocity 
hy 9 = ifc^f j[38j. Where the primes denote the derivatives with respect to conformal time rj, H is the conformal 
Hubble function and is the general sound speed, being defined as 

cl^'f. (11) 
bp 

From the intrinsic entropy perturbation T p2l . Wf\ , 



-r.(c^-cD^.^f^-^), (12) 

P VP P 



where the adiabatic speed of sound, c^, is purely determined by the EOS, 

2 P' 

c„ = — = w , 

" p' 'SH{l + w)' 



(13) 



it is found that only in the rest frame where S keeps gauge invariant can the general sound speed be a gauge invariant 
quantity since F and are gauge and scale independent. 

The following transformation gives the useful relation between the gauge invariant, rest frame density perturbation, 
(5, the density and velocity divergence perturbations in a general frame, S and 9 [22l. l37j: 

6 = S + 3nil + w)-^. (14) 

Thus we can express the pressure perturbation in a general frame, Sp, in terms of the general sound speed in the rest 
frame, c^, by using Eqs. (jl2ll4p 

Sp^clSp + 3m + w){cl-cl)p-^. (15) 
Using the above equation and Eq. (IT5|) . we have 

s' = -(1 + w){9 - 3(t)') - snicl - w)s - m{w' + m{i + w){cl - w))-^, (le) 
e' ^-n{i-3cl)9 + e f^^ + A . (17) 

\1 + w J 



^ Refer to the reviews of cosmological perturbation theory |36l44Clll for the detailed gauge transformation. 
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III. CONSTRAINT METHOD AND RESULTS FROM DATE SETS: WMAP, BAOS, AND SNIA 

A. Constraint Method and Data 

In the process of performing the observational constraints on cosmological parameters, we adopt the MCMC method 
by working with the publicly available CosmoMC package 41, 42], including the CAMB 43] code of calculating 
the theoretical CMB power spectrum. We modified the code for the RDE model with dark energy perturbations. In 
the first instance, work including perturbations of dark energy with EOS on either side of the cosmological constant 
boundary, i.e., w > — lorw < —1, was studied [i^]- In this case, the perturbed equations of dark energy can be solved 
well. However, when one encounters the perturbations of dark energy with EOS crossing the cosmological constant 
boundary -1, there inevitably exists a divergent problem due to the term 1 + w in the denominator in jj]q. (jl7[) . As 
pointed out in [l^, RDE behaves like a "quintom" when a < |. Here, we use the method proposed in |2^, dividing 
the whole range of w into three regions with a small positive parameter e: (1) w > — 1 + e; (2) — 1 + e > w > — 1 — e; 
and (3) w < —1 — e. Equations (|16ll7p are well satisfied in regions (1) and (3). The divergence happens in region (2). 
In order to deal with the divergence and keep dark energy perturbations continuous, the perturbations in region (2) 
are set as follows [2^ ]: 

S' = 0, 9' = 0. (18) 

While in the process of run ning the program, we set e — 10^^, which is within the limited range, e < 10~^[23]. By 
using the method, Refs. [2514281 ] gave the global fittings on the quintomlike dark energy with parametrized EOS when 
the dark energy perturbations were included. 

We adopt the following seven-dimensional cosmological spaces: 

P = {ui„LOa,es,T,a,ns,log[W'"A,]) (19) 

where ojf, = ili,h^ and Wc = flch^ are the physical baryon and cold dark matter densities, Qs is the ratio (multiplied 
by 100) of the sound horizon and angular diameter distance, t is the optical depth, a is the newly added Ricci 
parameter, rig is the scalar spectral index, and is defined as the amplitude of the initial power spectrum. The 
pivot scale of the initial scalar power spectrum we have used is fcso = 0.05Mpc~^. The corresponding priors are taken 
as follows: ^th"^ € [0.005,0.1], nch"^ € [0.01,0.99], Qs € [0.5,10], r e [0.01,0.8], a € [0.1,0.8], Hs € [0.5,1.5], and 
log[10^*'As] S [2.7,4]. In addition, we use a top-hat prior of the cosmic age, i.e., lOGyr < to < 20Gyr, impose a weak 
Gaussian prior on the physical baryon density flbh^ = 0.022 ± 0.002 from big bang nucleosynthesis [45|, and use the 
new value of the Hubble constant Hq = 74.2 ± 3.6km s~^Mpc~^ [i^ by a Gaussian likelihood function. 

In our calculations, we have taken the total likelihood L oc to be the products of the separate likelihoods of 

GMB, BAOs, and SNIa. Thus we have 

= XcMB + xIao + xlnia- (20) 

The full CMB data are used, including the new temperature and polarization power spectra from WMAP seven-year 
data [13] . We combine the additional distance measurements from BAOs and SNIa. For B AO information, we use the 
values of [r s{zd) / Dy -i) , rs{zd)/Dv{^-'i^)] and their inverse covariance matrix [1^. For SNIa, we use the 557 Union2 
SNIa data with systematic errors 5] . In order to assess the goodness of fit between the RDE model and ACDM model 
on the basis of the same observational data, we also perform the best fitting to the WMAP+BAO+SNIa data sets in 
the ACDM model. 



B. Fitting Results and Discussions 

In this section, we first present the constraint results from observational data. Then using the fitting results, we 
investigate the imprints on CMB temperature anisotropy, background evolution, and evolutions of perturbations for 
the RDE model with perturbations. 

The constraint results of full basic parameters and derived parameters are showed in Table HI where we list the 
mean values with Icr, 2a regions and best-fit values from WMAP alone in the RDE model, WMAP-f BAO+SNIa in 
the RDE model, and WMAP-|-BAO-|-SNIa in the ACDM model, respectively. It is seen that the inclusion of distance 
information from BAOs and SNIa leads to a significant change for best-fit value of parameter a. Both the best-fit 
values of the scalar spectral index, rig, in the RDE model are greater than 1. That is to say, the scalar spectrum 
is "blue" tilted. This result is different from that in the ACDM model, where within 2a regions the initial power 
spectrum with a blue tilt is not favored. For the combined constraint, the total best-fit of the RDE model is 
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Model and data 



RDE and WMAP alone 



RDE and WMAP+BAO+SNIa 



ACDM and WMAP+BAO+SNIa 



Parameters 



Mean 



Best fit 



Mean 



Best fit 



Mean 



Best fit 



Os 
r 
a 

iio 



log[10'°A, 



0.02221° 
0.12891^ 
1.0296t^ 
0.0986t^ 
0.2310t^ 
1.0159tJ! 

+0. 
0. 



3.113 



0006+0.0012 
0006-0.0011 
0064+0.0121 
0063-0.0125 
.0024+0.0047 
0023-0.0048 
0077+0.0309 
0092-0.0285 
0252+0.0526 
0257-0.0470 
.0186+0.0385 
.0186-0.0355 
036+0.073 
037-0.072 



0.0222 
0.1313 
1.0293 
0.0931 
0.2180 
1.0080 
3.109 



0.024lt° 
0.1087t;^ 
1.0302t° 
0.13741° 
0.34841° 
1.0887tg 

+0 




3.120 



0007+0.0014 
0007-0.0013 
0052+0.0104 
0055-0.0103 
.0024+0.0048 
.0025-0.0048 
0107+0.0417 
0122-0.0385 
0217+0.0433 
0215-0.0407 
.0181+0.0374 
.0185-0.0351 
048+0.100 
047-0.090 



0.0241 
0.1086 
1.0303 
0.1382 
0.3452 
1.0871 
3.122 



n n99c;+o ooo5+o.ooio 

U.UZZO_Q QQQ5_o 0010 
n 1 1 9/1 +0"033+0. 0065 
"•^^^^-0.0033-0.0063 
1 n'<Q9+00026+0.0051 
i.UOyz_o 0025-0.0050 



0.0877 



+0.0063+0.0245 
-0.0070-0.0230 



0.9680 



+0.0119+0.0238 
-0.0119-0.0235 
o 1 Q7+0. 036+0. 072 
J.lg)<_Q 037_o.071 



0.0226 
0.1123 
1.0393 
0.0877 

0.9691 
3.180 



fideO 

Age/Gyr 

fimO 

Ho 



0.7441^: 
14.00115; 

o.256j:E;- 

12.36li 
77.1313 



026+0.047 
025-0.056 
12+0.23 
11-0.22 
,025+0.056 
.026-0.047 
50+2.99 
50-2.96 
35+6.56 
31-6.57 



0.745 

14.03 
0.255 
12.02 
77.54 



13.83tg; 
0.259t^ 
14.49ti- 
71.63t} 



017+0.031 

017-0.034 

12+0.23 

12-0.23 

017+0.035 

017-0.031 

68+3.28 

69-3.41 

66+3.30 

63-3.21 



0.746 
13.82 
0.254 
14.61 
72.26 



p. r,9r,+ 0. 014 + 0. 027 

' -0.015-0.03 
1 •? 7C+0. 011+0. 021 
iO. <0_o. 011-0. 021 
r, 970+0.015+0.03 
U.z/O_o 014-0.027 



10.51 

70.32 



+ 1.18+2.34 
-1.17-2.35 
+ 1.28+2.52 
-1.27-2.47 



0.728 
13.75 
0.272 
10.32 
70.38 



TABLE I: The mean values with la, 2a regions and best-fit values for the RDE model from WMAP alone and the combination 
of WMAP, BAOs, and SNIa, respectively, and the mean values with la, 2a regions and best-fit values for the ACDM model 
from the combination of WMAP, BAOs and SNIa. 



8056.1, while = 8002.6 for the best-fit ACDM model. Therefore, the ACDM model still shows a better fit to 
the current data than the RDE model. In Fig. [TJ we give the one-dimensional (ID) marginalized distributions of 
parameters and 2D contours with confidence level (C.L.) from the combination of WMAP, BAOs, and SNIa. The 
respective contributions of WMAP, BAOs, and SNIa to cosmological parameters can be seen in Fig. [5J which shows 
the 2D contours of parameters (a, Vl,no) from WMAP alone, BAOs alone, SNIa alone, and their combinations. The 
stringent results can be obtained from the joint constraints of dynamical and geometrical aspects. 

In Fig. 131 we show CMB temperature power spectra for the RDE model by using the best-fitting results in Table |T] 
and for the ACDM model with the best-fit values from the WMAP alone constraint [H and WMAP+BAO+SNIa in 
Table HI It is found that there are significant discrepancies of CMB temperature power spectra in the two models on 
large angular scales, where the contribution to the secondary CMB temperature power spectrum is mainly from the 
late ISW effect. Both the curves of the CMB temperature power spectra for the RDE model on large angular scales 
almost tend to be plat. A large-scale [l < 10) plateau happens in the RDE model. In order to understand the imprint 
on the CMB temperature power spectrum from parameter a alone, we illustrate the CMB temperature power spectra 
for different a in Fig. H) It is shown that the CMB temperature power spectra on all the scales decrease as the value 
of a becomes smaller. Because a has an effect on the expansion rate, the angular diameter distance to recombination 
becomes larger when a decreases. So there is the rightward shift of the positions of peaks. What is more, the usual 
dimensionless densities of the RDE and matter component are changed for different parameters, a, as shown in Fig. 
[SJ That is to say, the change in the expansion rate brings about the differences of usual dimensionless densities when 
there is the same matter density, ^mo- We can see at early times the usual dimensionless matter density decreases 
when a become larger. We try to look into the role of the change of a by taking another point of view. We consider 
that this effect can be viewed as coming from the different effective matter densities based on the same expansion 
rate, i.e., + zY /El. So there is a lower effective matter density with the increase in a. The lower effective 

matter density causes the delay of the epoch of matter-radiation equality, i.e., the smaller z^q. This change in redshift 
of equality leads to the enhanced late ISW effect on large scales and power spectra on small scales. Furthermore, 
the delay of the time of dark energy domination for smaller a also leads to the suppressed late ISW effect on large 
scales. In Fig. [51 we show the evolutions of cold dark matter and dark energy perturbations on the large scale, i.e., 
k = 0.001Mpc"\ 

The differences of background evolutions between the ACDM model and RDE model can be seen through the 
evolutions of the distance moduli and relative distance moduli in Fig. [T] Then following the method of the propagation 
of errors by the Fisher matrix analysis [49l - l5T| ^, we plot the evolution of EOS with Icr errors in Fig. |8]by using the 
mean values from WMAP+BAO+SNIa. The current value of the EOS is -1.0853, less than -1. What is more, it is 



^ We refer to [s^ to marginalize the cosmological parameters which are not associated with the EOS and obtain the error covariance 
matrix of EOS parameters from the sub-Fisher matrix. 
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FIG. 1: The ID marginalized distributions on individual parameters and 2D contours with 68%C.L., 95 %C.L., and 99.7%C.L. 
between each other using the combination of the observational data from WMAP, BAOs and SNIa for the RDE model. 



seen that the EOS of RDE crosses the cosmological constant boundary at z ^ 0.05. 



IV. SUMMARY 



In summary, we perform a global fitting on the RDE model with perturbations to the combined constraints from 
the full CMB data, BAOs, and SNIa by using the MCMC method. All the cosmological parameters for the RDE 
model are well determined, as shown in Table HI It is worth noting that the initial power spectrum with a blue tilt is 
favored for the RDE model. Making use of the constraint results, we investigate the CMB temperature power spectra. 
It is found that the differences of CMB temperature power spectra on large angular scales are significant between the 
ACDM model and the RDE model with perturbations. The main feature of CMB temperature power spectra on large 
angular scales for the RDE model with perturbations is that there exists a flat plateau. In view of the differences of 
large-scale power spectra, on which the ISW effect plays an important role, the effects on RDE parameters made by 
the ISW data from the cross correlations of CMB andlarge-scale structure will be worth investigating in future work. 
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0.35 




0.15 - 



0.15 



0.25 



0.35 



0.45 



0.55 
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FIG. 2: The 2D contours with 68%C.L., 95 %C.L., and 99.7%C.L. between a and Q.^o for the RDE model from WMAP alone 
(red dashed line), BAOs alone (magenta dash-dotted line), SNIa alone (blue dotted line), and the combination of WMAP, 
BAOs and SNIa (solid black line). 
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FIG. 3: CMB temperature power spectra vs multipole moment I, where the black dots with black error bars denote the 
observed data with their corresponding uncertainties from WMAP, the red solid line is for the ACDM model with the best-fit 
values from the WMAP alone constraint Ij, the green dash-dotted line is for the ACDM model with the best- fit values from 
the WMAP-f BAO-f SNIa joint constraint, the blue dashed line is for RDE with the best-fit values from the WMAP alone 
constraint, and the magenta dotted line is for the RDE model with the best-fit values from the WMAP+BAO+SNIa joint 
constraint. 




FIG. 4: CMB temperature power spectra vs multipole moment I for the RDE model with different parameters, a, where the 
blue dotted line is the power spectrum with a — 0.44, the black solid line is the power spectrum with the best-fit value of a 
from the WMAP+BAO+SNIa joint constraint, and the magenta dashed line is the power spectrum with a — 0.24. 
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FIG. 5: The usual dimensionless energy densities of RDE (thin) a-nd matter component (thick) ^ ^"^o(^+^) ^ ]-e(jg}jif); 

z for the RDE model with different parameters, a, where the blue dotted lines are dimensionless energy densities with a — 0.44, 
the black solid lines are dimensionless energy densities with the best-fit value of a from the WMAP+BAO+SNla joint constraint, 
and the magenta dashed lines are dimensionless energy densities with a = 0.24. 




r| / Mpc 



FIG. 6: The evolutions of cold dark matter density perturbation, 5c (dash-dotted line), dark energy density perturbation, \5de\ 
(solid line), and its velocity perturbation, Vde (dashed line). 
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